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We investigate the quantum hadrodynamic equation of state for neutron stars (with and without 
including hyperons) in the presence of strong magnetic fields. The deduced masses and radii are 
consistent with recent observations of high mass neutron stars even in the case of hyperonic nuclei 
for sufficiently strong magnetic fields. The calculated adiabatic index and the moments of inertia 
for magnetized neutron stars exhibit rapid changes with density. This may provide some insight 
into the mechanism of star-quakes and flares in magnetars. 

PACS numbers: 



I. INTRODUCTION 

Soft 7-ray repeaters (SGRs) and anomalous X-ray 
pulsars (AXPs) are believed to be evidence for magne- 
tars, i.e. neutron stars with surface magnetic fields of 
10 14 ~ 10 15 G In the interior of these magnetic neu- 
tron stars, the magnetic field strength could be as much 
as 10 18 G according to the scalar virial theorem. Such 
strong magnetic fields may affect the properties of neu- 
tron stars such as the relative populations of various par- 
ticles, the equation of state (EOS), and the mass-radius 
relation. Many studies of dense nuclear matter in the 
presence of strong magnetic fields have been reported @- 
9]. These works have considered the electromagnetic in- 
teraction, the Landau quantization of charged particles, 
and the baryon anomalous magnetic moments (AMMs). 
However, a detailed analysis of neutron stars with strong 
magnetic fields is still an area of active research. 

Recently [l(il |. a giant 7-ray flare, SGR 1806-20, has 
been observed. The total flare energy was estimated to be 
as much as 2 x 10 46 erg. This is ~ 10 2 times higher than 
the two previously observed giant flares [ill El > an d is 
believed to have been released during a reconfiguration of 
magnetic fields of the neutron star. In Ref. |13| . the sim- 
ilarities between SGR events and star-quakes relating to 
the sudden change of pulsar periods were discussed, and 
the possibility that SGRs may be powered by starquakes 
was suggested. These phenomena might be explained by 
sudden changed in the magnetic pressure. 

Although many uncertainties remain regarding star- 
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quakes and the strength of magnetic fields in the interior 
of neutron stars, it has been suggested [|| that neutron 
star matter might be unstable in the presence of strong 
magnetic fields because of the onset of discrete Landau 
level energy quantization. This can cause rapid changes 
in the pressure response to changes in density. This in- 
stability could be one source for star-quakes. Thus, one 
may conjecture the following unifying scenario from star- 
quakes to pulsar glitches. If a strong magnetic field can 
cause star-quakes, it may crack the surface of the star. 
Magnetic field energy may then be released through the 
cracks. The released magnetic energy and associated re- 
connection may be observed as a SGR or a giant flare. In 
addition, the release of the magnetic field energy density 
may affect the equation of state (EOS) of the neutron 
star. This, in turn, could give rise to a change of the mo- 
ment of inertia, and thus, the neutron-star spin period 

In Ref. [ij properties of magnetic stars were studied, 
but only in the context of a cold free n,p,e Fermi gas. 
It could not be determined in that work whether the 
features of interest would persist with a realistic nuclear 
equation of state. It is useful, therefore, to reconsider the 
structure and dynamics of magnetic neutron stars in the 
context of a realistic nuclear equation of state. Therefore, 
in this paper we calculate the populations of particles, the 
instability in the adiabatic index F, the EOS, and the 
moment of inertia for various magnetic field strengths by 
using the techniques of quantum hadrodynamics (QHD). 
We then discuss the possibility that star-quakes and SGR 
flares might be explained by the rapid change of the adi- 
abatic index T with density due to the population of 
Landau levels. We find that, similar to the Fermi gas 
approximation Q, the magnetic QHD model also shows 



2 



an instability, i. e. a sudden change of the adiabatic index 
as the density increases. As in the free nucleon gas this 
is attributable to the discrete excitation of Landau levels 
as the density increases, but also to the appearance of 
hadronic species at high density. These changes in the 
adiabatic index lead to a sudden change in the pressure 
response of the star. Therefore, star-quakes and an as- 
sociated release of magnetic field energy may take place. 
Moreover, when we assume that magnetic fields as large 
as 10 18 G exist inside the star, we find that a release of 
magnetic field energy could decrease of the moment of 
inertia leading to an increase in the the spin rate of the 
star. 

In section II, we briefly introduce our theoretical 
framework for magnetized dense matter based upon the 
QHD approach. A method for calculating the possible 
change in the moment of inertia by the abrupt variations 
of the adiabatic index is also presented. Discussions of 
star-quakes and the neutron-star spin are presented along 
with numerical results in section III. A summary and con- 
clusions are given in section IV. 



II. THEORY 

A. Relativistic mean field Lagrangian with strong 
magnetic fields 

The Lagrangian of the QHD model for dense nuclear 
matter in the presence of a magnetic field can derived 
by introducing of a vector potential . The resulting 
Lagrangian can be written in terms of the baryon octet, 
leptons, and five meson fields as follows 
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where the indices b and I denote the baryon octet and the 
leptons (e~ and p~), respectively. The effective baryon 
mass, M & *, is simply given by M 6 * = M b - g ab a - g a *bcr*, 
where Mb is the free mass of a baryon in vacuum and the 
g a b are associated coupling constants. The <r, w and p 
meson fields describe the nucleon-nucleon (N — N) and 
nucleon-hyperon (N — Y) interactions. The Y — Y inter- 
action is mediated by the a* and <fi meson fields. U(a) 



is the self interaction of the a field given by U(a) = 
\92^ + j93<7 4 ■ W^, Ri pv , and F^ v denote the field 
tensors for the w, p, <f) and photon fields, respectively. 

The anomalous magnetic moments (AMMs) of the 
baryons interact with the external magnetic field via 
terms of the form of Kb<J^ v F^ 1 ' where = §[7 M ,7„] 
and Kb is the strength of the AMM for a baryon, i.e. 
K p = 1.7928/iAr for a proton with p^ the nuclear magne- 
ton. 

In general, the AMMs could depend upon the matter 
density. Therefore, one can take account of the medium 
effect through density dependent AMMs which can be 
evaluated within the quark meson coupling (QMC) model 
[9j. In this report, however, we did not take account 
of the effect of density dependent AMMs because the 
calculation here is performed within the QHD model. 

Energy spectra for the baryons and leptons are given 

by 
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where E^ and E^ denote the energies of charged and 
neutral baryons, respectively. The Landau quantization 
for charged particles in a magnetic field is denoted as 
^ = n+ l/2 — sgn(q)s/2 = 0, 1, 2 • • • , where the sign of 
the electric charge is denoted as q sgn(q) and s = 1(— 1) 
is for spin up (down). 

Chemical potentials for the baryons and leptons are 
given by 



p b = E f + g^b^o + 50600 + g p b! b z P3o, 
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where Ej is the baryon Fermi energy and fey is the lepton 
Fermi momentum. For charged particles, £y is written 

as 



Ef = kf 



m*^ + 2v\q b \B 
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for 



where fc* is the baryon Fermi momentum and v 
neutral baryons. 

We apply three constraints for calculating the proper- 
ties of a neutron stars: 1) baryon number conservation; 
2) charge neutrality; and 3) chemical equilibrium. The 
meson field equations are solved along with the chemical 
potentials for baryons and leptons subject to these three 
constraints. The total energy density is then given by 
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£tot — £m + st, where the energy density for the matter 
fields is given as 

b I 

+ \™lu 2 + \™\<i> 2 + + U(a), (6) 

and the energy density due to the magnetic field is given 
by ef — B 2 /2. The total pressure is then 



Ptat = P m + ~B Z 



(7) 



where the pressure due to the matter fields is obtained 
from P m — ^if-iiPl, — £ m - The relation between the 
mass and radius for a static and spherical symmetric 
neutron star is generated from a solution to the Tolman- 
Oppenheimer-Volkoff (TOV) equations using the equa- 
tion of state (EoS) described above. 



B. Slowly rotating neutron stars 

To calculate the moment of inertia for a slowly rotat- 
ing neutron star, we follow the methods detailed in Refs. 
14 , [Til l which are briefly summarized here. The metric 
of an axially symmetric neutron star can be written in 
generalized rotating Schwarzschild coordinates as 



ds 2 
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If the neutron star is rotating uniformly with a stellar 
frequency Q far below the Kepler frequency 



(9) 



the slow-motion approximation is valid and the moment 
of inertia can be written: 



J 
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where J is the angular momentum, while v{r) and w(r) 
are radially-dependent metric functions. M(r), and 
P(r) are the mass of the star, energy density, and pres- 
sure, respectively, derived from a solution to the TOV 
equation. 

The metric functions v{r) and A(r) in Eq. © are 
unchanged from the values for a spherically symmetric 
neutron star. Thus, A(r) is obtained from the mass M(r) 
by the usual Schwarzschild solution: 



<?ii(r) = e 2A « = (l-2M(r)/r)' 
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and v{r) can be calculated from 
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The metric function w(r) denotes a frequency corre- 
sponding to the dragging of local inertial frames by the 
rotating star. The relative frequency Co(r) = 51 — uj(r) 
appearing in Eq. (fTU|) represents the angular velocity 
of the fluid as measured in a local inertial reference 
frame. In particular, the dimensionless relative frequency 
w(r)=w(r)/f2 satisfies the following second-order differ- 
ential equation 
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where 



j(r) = e 



y/l - 2M(r)/r if r < R , 
1 if r > R . 

(14) 



Note that to(r) is subject to the following two boundary 
conditions: 
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In numerical calculations, one can integrate Eq. (|13[) from 
the center of the star to its surface for various central fre- 
quencies. When the boundary conditions at the surface, 
Eq. (|16p are not satisfied for an arbitrary choice of u5 c , 
one must rescale the function and its derivative by an 
appropriate constant to correct for the mismatch. After 
solving both u)(r) and /, one can check the consistency of 
the results by testing the accuracy to which the following 
condition is satisfied: 



u?(R) 



61 

R 4 



(17) 



III. RESULTS AND DISCUSSION 



We use the parameter set given in Ref. [l6[ for the 
coupling constants, ga-N, g^N and g p N, where N denotes 
the nucleon. The coupling constants for hyperons in a nu- 
clear medium, g u y and g G y , are determined by the quark 
counting rule and the relevant hyperon potentials at the 
saturation density. The strength of these potentials is 
fixed at U A = -30 MeV, J7 S = 30 MeV and U s = -15 
MeV. Since the magnetic fields may also depend upon 
the density, we utilize the form for the density-dependent 
magnetic fields suggested in Refs. @, H, Q 

B (p/p ) = B™* + B [1 - exp{-/3 (p/po) 1 }] , (18) 
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where B sur f is the magnetic field at the neutron-star 
surface, taken from observations to be ~ 10 15 G, and J5 
represents the magnetic field saturation strength in the 
high density region. 

In this work, we use /3 = 0.02 and 7 = 3. Since the 
magnetic field is usually written in units of the critical 
field for the electron B c e = m 2 e /e = 4.414 x 10 13 G, the 
B and the B in Eq. (|T8|) can be expressed as dimen- 
sionless quantities B* = BjB c e and Bq = Bq/B^. Based 
upon this magnetic field representation, we investigate 
the structure of neutron stars as a function of magnetic 
field strength both with and without hyperonic matter. 
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A. Particle populations, adiabatic index and 
equation of state 

The populations of particles as a function of density 
in our model are shown in Fig. [TJ for four cases. Shown 
are particle fractions Yi for both pure nucleonic and hy- 
peronic equations of state, and both with and without 
magnetic fields. Since the details about the effects of 
magnetic fields in neutron stars have been discussed in 
our previous paper Q, we briefly explain the magnetic 
field effects on populations of particles. In Fig. [U one 
can notice the effects of a magnetic field from the differ- 
ence between the (no B) and (Bq = 5 x f 4 ) figures for 
both nucleonic (upper) and hyperonic phases (lower). 

In general, charged particles are strongly affected by 
the EM interaction term (eB), which gives rise to an 
increase in the electron fraction. Therefore the condi- 
tion of charge neutrality and the EM interaction enhances 
the proton fraction. This suppresses hyperons by baryon 
number conservation. This phenomenon explicitly ap- 
pears if we utilize Bq = 10 5 as in our previous paper [9(. 
But its effect becomes indiscernible in the present case 
with a weaker field strength, Bq = 5 x 10 4 . 

The most prominent feature caused by the strong mag- 
netic field in Fig. 1 is the kink pattern in the populations 
of electrons and muons, which is clearly shown for the 
hyperonic phase. Of course, that pattern is caused by 
the change of Landau levels owing to the existence of a 
magnetic field. 

In Fig. [21 we show the adiabatic index V = = 

(1 + f>)4j" f° r both the nucleonic (upper) and the hy- 
peronic (lower) equation of state. This is of particular 
interest because the adiabatic index may be used as a 
critical factor for understanding the radial stability of a 
star [13, El. 

In the nucleonic phase, one can see multiple kinks sim- 
ilar to the results shown in Ref. [3| based upon a cold 
Fermi gas model. This oscillatory pattern takes place 
as charged particles begin to rapidly fill the next unoc- 
cupied Landau level when the density increases. In the 
hyperonic phase, one can see a sudden large increase in 
the adiabatic index at ~ 3po- This is caused by the pro- 
duction of S~ hyperons around that density as shown in 
Fig. I. 
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FIG. 1: Populations of particles in a neutron star. The upper 
two panels are for the nucleonic phase and the lower two pan- 
els are for the hyperonic phase with and without a magnetic 
field, respectively, where we used Bq — 5 x f 4 as explained 
in the text. Note, this field strength is a bit different value 
from the value used in our previous paper [9(. 
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FIG. 2: Adiabatic index for the nucleonic (upper) and hyper- 
onic (lower) equations of state. The adiabatic index is defined 
asr= = (1 + Jr)^. Here Bq is given in units of 10 4 . 
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FIG. 3: Equation of state pressure as a function of density in 
units of the nuclear saturation density. 



In the high density region beyond 3/?o, similar kink pat- 
terns appear in both the nucleonic and hyperonic phases. 
At high density, however, the kinks are relatively smaller 
than the kinks in the lower density region. This is be- 
cause the high density increases the occupied levels but 
does not cause a rapid shift in the Landau levels. This 
means that there should exist a critical density for such 



discontinuous pressure response, at least in the hyperonic 
phase. 

Similarly to the Fermi gas model Q, the magnetic 
fields in the QHD model also give rise to sudden jumps 
from the occupied Landau levels to the unoccupied lev- 
els. This leads to kinks in the adiabatic index because of 
the unstable structure in the magnetic pressure. These 
oscillatory kinks can be conjectured to cause star-quakes. 
Specifically, in the hyperonic phase, the instability caused 
by the magnetic field becomes significantly larger than 
that in the nucleonic phase. 

In Fig. [3l the EOSs are shown for both nucleonic and 
hyperonic matter. Magnetic fields can make the EOS 
stiff in both phases. In particular, since the magnetic 
fields suppress hyperons, the effect of stiffness owing to 
the magnetic field in the hyperonic phase is far larger 
than in the nucleonic phase. Thus, the stiff EOS gives a 
large mass and radius for a neutron star. That is shown 
in Fig. |U 



B. Mass-radius relation and the moment of inertia 

The mass-radius relations for both the nucleonic and 
the hyperonic phases are compared with observational 
limits in Fig. 4. In this context, it is particularly note- 
worthy that larger maximum masses and radii are possi- 
ble when a strong interior magnetic field is present. Neu- 
tron stars and heavy ion collision data provide valuable 
constraints on the EOS for dense matter There are 
indeed recent reports of evidence for higher maximum 
masses for neutron stars. For instance, M — 2.O±O.1M0 
for 4U 1636-536 was reported in Ref. [20}. Very recently 
authors in Ref. [2l[ investigated seven neutron stars, six 
binaries and an isolated neutron star (RX J1865-3754) 
and showed that M « 1.9 - 2.3M© and R « 11 - 13 km. 
That is labeled as the SLB data in Fig. 4, which also 
shows the 2a lower and upper limits. 

Pulsar I of the globular cluster Terzan 5 (Ter 5 I) shows 
a lower mass limit Mm. star > 1.68Af Q at the 95 % con- 
fidence level (22j |. Another constraint deduced indepen- 
dently of given models is obtained from XTE J1739-285 
[23[ | . That work provides a constrained curve for the ratio 
between mass and radius. Although we do not indicate it 
here, another recent paper [24| has also reported a mass 
of M « 1.97M© for PSR J1614-2230 from a detection of 
the Shapiro time delay. 

For the case of a nucleonic EOS, the mass and radius 
relations satisfy all constraints. The hyperonic phase, 
however, shows some discrepancies. Although a hyper- 
onic phase cannot satisfy all of the constraints, magne- 
tized hyperonic neutron stars can explain large neutron- 
star masses within the constraints from the observational 
data. 

In Fig. [5j we show the moment of inertia in units of 
J45 = 10 45 g cm 2 in order to investigate the spin up of 
the star. One can easily notice the large difference in 
the moment of inertia between stars with and without 
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FIG. 4: Mass-radius relations for nucleonic (upper) and hy- 
peronic (lower) equations of state. Three constraints are 
shown in the figure. XTE is a short expression of XTE J1739- 
285 [23| , Ter 5 I is Pulsar I of the globular cluster Terzan 
5(M st ar = 1.68M Q ) [22], and SLB is taken from the data 
(2cr limits) on r p h S> R in Ref. [2l|] where r p u is the photo- 
spheric radius. values are given in units of 10 4 . Detailed 
explanations are given in text. 



magnetic fields. This implies that reducing the magnetic 
field strength could cause a decrease in the moment of 
inertia. If there is a mechanism for this reduction, one 
might expect an decrease in the spin period by angular 
momentum conservation. 

For illustration, we plot ^45 as a function of Bq on Fig. 
[5] for models with a central density of p = 4/?o- For a can 
field strength around Bq = 10 4 , relatively small changes 
of the magnetic field give rise to a significant change in 
the moment of inertia leading to a spin change of the 
neutron stars. 

To examine the effect of SGR flares on the moment of 
inertia, consider the released magnetic energy is Efi are ~ 
2 x 10 46 erg for the largest observed magnetar flare yet 
observed [lol ]. If we assume that this magnetic energy 
released from a region below the surface of thickness of 
thickness Ar, than the released energy is related to the 
change in the magnetic energy density (Ai3) 2 /2) times a 
volume factor, Efi are = (47ri? 2 Ar) x (AB) 2 /2 where R 
is the radius of the emission region taken to be ~ 10 km. 
An emission ~ 10 46 erg corresponds to a change in the 
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FIG. 5: The moment of inertia (745 = 10 45 g • cm 2 ) for nu- 
cleonic and hyperonic equations of state with and without 
magnetic fields as a function of of radius (upper) and mag- 
netic field strength (lower) figure. 



magnetic field of only AB ~ 10 16 G from a thickness of 
Ar = 1 m ( or - 10 15 G for Ar = 100 m ). From Fig. [5] 
it is apparent that a change of ~ 10 16 G corresponds 
to about a 1% decrease in the moment of inertia and 
the rotation period. For a typical flare associated with 
AXPs, however, the energy released is only ~ 10 40 erg 
[25l ] . corresponding to AB ~ 10 13 G, implying a decrease 
of the rotation period by a factor of ~ 10 -5 . Interest- 
ingly, this is comparable to the glitches observed f2(| in 
some AXPs in association with a radiative event. Thus, 
an association of AXP emission and glitches is possible 
in this mechanism. 



IV. SUMMARY 

We have investigated neutron stars which contain 
strong interior magnetic fields. The populations of par- 
ticles, the adiabatic index, the EOS, the mass-radius re- 
lation, and the moment of inertia have been calculated 
by using the QHD Lagrangian modified to include strong 
magnetic fields. We discuss the effects of a strong interior 
magnetic field on the physical quantities characterizing 
neutron stars or magnetars. 
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Possible changes in the mass-radius relation for neu- 
tron stars with strong magnetic fields are compared with 
the constraints deduced from observational data. The 
mass-radius relations for both nucleonic and hyperonic 
phases with strong magnetic fields satisfy the maximum 
mass suggested by Ter 5 I and XTE J1739-285, but the 
constraint in Ref. [2lj turns out to be inconsistent with 
the results for the hyperonic phase unless a strong inte- 
rior magnetic field is present. 

The possibility of star-quakes, magnetar flares and 
changes in the spin rate of AXPs associated with flares 
has also been addressed by determining the adiabatic in- 
dex and the moment of inertia with a magnetic QHD 
equation of state. Our results suggest that changes in 
the pressure response of the star with changes in density 
could lead to star quakes and magnetic flare energy re- 
lease consistent with magnetar flares if the interior mag- 



netic field strength is large enough. The implied change 
in the moment of inertia associated with magnetic flares 
is also consistent with glitches in the AXP period occa- 
sionally associated with radiative events. 
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